Emergent Gravity at a Lifshitz Point from a Bose Liquid on the Lattice 
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We propose a model with quantum bosons on the fee lattice, which has a stable algebraic Bose 
liquid phase at low energy. We show that this phase is described by emergent quantum gravity at 
the Gaussian z — 3 Lifshitz fixed point in 3 + 1 dimensions. The stability of this algebraic Bose 
liquid phase is guaranteed by the gauge symmetry of gravitons and self-duality of the low energy 
field theory. By tuning one parameter in the lattice boson model we can drive a phase transition 
between the 2 = 3 Lifshitz gravity and another algebraic Bose liquid phase, described by gravity at 
the z = 2 Lifshitz point. 



The challenge of finding a satisfactory theory of quan- 
tum gravity has stimulated theoretical physics for many 
decades. We expect that such a theory should make sense 
of the quantum fluctuations of the spacetime metric at 
low-energy scales, while providing a quantum mechanical 
completion of the system at short distances. Recently, a 
new approach to this long-standing puzzle has been pro- 
posed In this approach, gravity is treated using the 
traditional path integral methods of quantum field the- 
ory, but without assuming Lorentz invariance as a funda- 
mental symmetry at short distances. The gauge symme- 
tries are those spacetime difFeomorphisms that preserve 
a preferred foliation of spacetime by fixed time slices, 
generated by 



6x'=C{x^,t), St = C{t). 



(1) 



In terms of the spatial metric gij{x^,t), the shift vector 
Ni{x^ ,t) and the lapse function N{t), the action is 

Lk - 5]A'.,i^^^ -A(^i^:)2, (2) 



where 



D,N, - D,N,) 



(3) 



is the extrinsic curvature of the constant time slices in the 
spacetime foliation, Di is the covariant derivative defined 
by gij, K and A are coupling constants, and the Ly can be 
an arbitrary local Lagrangian built from gij , its Riemann 
tensor jke, and the covariant derivative Di, without the 
use of time derivatives. 

In this broader framework, new Gaussian fixed points 
of the renormalization group (RG) are possible. These 
novel fixed points are characterized by a dynamical scal- 
ing exponent z ^ 1 (which depends on the choices made 
in Lv), and they exhibit a Lifshitz-type scaling. Given 



the possibility of such Lifshitz gravity fixed points, it is 
natural to ask whether they can serve to provide an ul- 
traviolet (UV) completion of gravity, and whether the 
pattern of the RG flow can restore z = 1 and Lorentz 
invariance in the infrared (IR) regime, where the theory 
can be tested against general relativity. 

In this paper, we show that this new class of Lifshitz 
gravity fixed points P, Q with z — 2 and z = 3 in 3 + 1 di- 
mensions can also emerge as infrared fixed points, in the 
low-energy limit describing certain lattice systems. The 
microscopic degrees of freedom of these systems are of 
the conventional type familiar in condensed matter, but 
their collective behavior leads at low energies to novel 
gapless phases, described by the Lifshitz gravity fixed 
points. The microscopic lattice degrees of freedom ex- 
hibit no gauge symmetry: The foliation-preserving gauge 
invariance at the Lifhitz gravity fixed points is entirely 
emergent in the low-energy description of the lattice sys- 
tem. The more fundamental physical objects on short 
length scales give rise to the Lifshitz gravitons as their 
low-energy and long-distance collective excitations. 

An analogous mechanism, with an emergent U(l) 
gauge symmetry, was realized previously in the quan- 
tum dimer model on the cubic lattice. In this model, the 
short distance excitations are spin singlet valence bond 
fluctuations, while the long length-scale excitations are 
photons [1, 01 ■ A similar photon phase can also be con- 
structed in spin models on the pyrochlore lattice 0, § . 
Given these examples, it is natural to look for a lattice 
model whose infrared behavior is controlled by the Lif- 
shitz gravity fixed points of [1, 2]. 

It should be noted that there are two distinct ways in 
which one can attempt to relate Lifshitz gravity to lattice 
models. In this paper, we work on a fixed rigid lattice, 
and find degrees of freedom whose long distance dynam- 
ics is captured by the Gaussian fixed points of Lifshitz 
gravity. The lattice is nondynamical. Alternatively, one 
can try to obtain Lifshitz gravity as a continuum limit of 
a lattice model defined as a sum over a suitable class of 
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random triangulations of spacetime geometries. Here it is 
the lattice itself that is dynamical, and no extra degrees 
of freedom are invoked. The most promising candidate 
for such a random lattice model is offered by the causal 
dynamical triangulations (CDT) approach to quantum 
gravity. In the CDT approach (see Ref. for a review), a 
summation over random lattices constrained to respect a 
preferred foliation structure of spacetime serves as a non- 
perturbative definition of quantum gravity, and yields a 
continuum limit with four macroscopic spacetime dimen- 
sions at long distances. It has been suggested in i (see 
also the recent paper that the CDT approach might 
be viewed as a lattice regularization of Lifshitz gravity. 
Further evidence for this scenario comes from the qual- 
itative behavior of the spectral dimension of spacetime, 
which indicates that the model flows from a z = 3 UV 
fixed point to an z = 1 fixed point at long distances 0]. 

In condensed matter systems, gapless bosonic exci- 
tations are usually Goldstone modes of certain sponta- 
neously broken continuous symmetry. For instance, the 
phonon modes of solids are the Goldstone modes of trans- 
lation symmetry, and the magnons are the Goldstone 
modes of the broken spin rotation symmetry in the mag- 
netic ordered phase. There has been a lot of interest 
recently in finding bosonic phases with gapless excita- 
tions which do not originate from breaking a symmetry, 
and which do not require fine tuning. Such phases are 
referred to as the algebraic Bose liquid (ABL) phases, be- 
cause the boson density correlation falls off algebraically 
in space and time. An ABL phase is defined and char- 
acterized by its low-energy field theory. Searching for 
stable ABL phases, or more generally Bose metal phases, 
is one of the active fields in condensed matter theory 
The photon phase of the quantum dimer model 
mentioned above is a well-understood stable ABL phase 
which was first studied in high energy community and 
then widely applied to condensed matter systems, espe- 
cially fractionalized states of strongly correlated systems 

SO- 

Can we find other examples of ABL phases, especially 
in three spatial dimensions? In Ref. 13, Hi, a novel ABL 
phase with z — 2 dispersion was proposed in a quan- 
tum boson model on the face centered cubic (fee) lattice 
with only local boson hopping and density repulsion. Al- 
though the microscopic model only has the lattice point 
group symmetry and the U(l) global symmetry corre- 
sponding to the conservation of the total boson number, 
a gauge symmetry similar to linearized diffeomorphisms 
emerges at low energies. We will review the low-energy 
effective field theory of this model below, and show that 
it is given by the Lifshitz gravity of Ref. [ij at the z = 2 
Gaussian fixed point. In the lattice construction, the 
z — 2 dispersion is protected by the gauge symmetry 
which emerges at low energy and by the microscopic dis- 
crete symmetries of the underlying degrees of freedom, 
and the stability of this ABL phase is guaranteed by its 



self-duality. In this work, we show that by turning on 
one extra density repulsion term in the model proposed 
m Ref. ^E^, one can drive a phase transition between 
the z = 2 phase mentioned above and another stable 
z — 3 phase. We will show that this z — 3 phase is also 
a stable ABL phase, with a self-dual structure at low en- 
ergy. The field theory of this ABL phase is identical to 
the z — 3 Lifshitz gravity of Ref. [2] at the Gaussian fixed 
point. 

We start with describing the full Hamiltonian of our 
lattice boson model. This model is defined on the fee 
lattice: The physical quantities will be defined on both 
the sites and the unit square faces of a cubic lattice. We 
denote each site of the cubic lattice by r = (r^ ,ry,rz), 
and each unit square in the plane by f ± | ± |, 

with i,j = x,y,z. As our dynamical variables, we as- 



sign three boson numbers {nxx,ny 



2) on each site of 



a cubic lattice, and one boson number riij to each face in 
ij plane of the cubic lattice. The corresponding creation 
and destruction operators wiU be denoted by bl^,bii,blj 
and bij. The microscopic Hamiltonian contains the fol- 
lowing terms, 

H 



Ht 



Hf + H„ + H,, + H„ 



-tiHsp — t2Hp 



X link — V {2nxx,f ~^ '^^xx.f+x ^ 



xy.r-\- 1 



"xZjT+^ + i 



''xz,r+? 



8n)2, 



EU2 
2 

i<j 



Ht include all the local boson hoppings. H^p is the 
hopping between each site and its adjacent plaquettes, 
Hsp = 4fe,Aj,r±|±l while Hpp is the hop- 

ping between two nearest neighbor plaquettes (Fig. [T]i). 
The exact values of the amplitudes ti , ^2 are unimpor- 
tant, and we will tentatively take both of them at the 
same order of magnitude ti , t2 ~ t. Hy is a large density- 
density interaction between bosons, Hy^ ^ link in Eq- © 
is the part of Hy associated with one link along the x 
direction; it is a sum of the boson numbers on the two 
sites and the four faces sharing this link, and n is the 
average boson filling on each quantum state. Contribu- 
tions to Hy for links along y and z directions are defined 
analogously. H^ is a small repulsive interaction between 
each fiavor of bosons on both sites and faces, with ui 
generally not equal to U2. Finally, Hyi is another on- 
site repulsive interaction on each site of the cubic lattice. 
The lattice structure and the distribution of the boson 
number is shown in Fig. [TJi. 

If we first take V = 0, this model will reduce to the 
model constructed in Ref. n Hy is dominant. 
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the components of Aij transform as 
T : i y t, Aij y -^ijj 
Tf. ■ f )■ r + fc, Aij >■ Aij^ 
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FIG. 1: (a), The lattice structure of this model, with the 
coordinates of site, link and plaquette. The term j unk in 
Eq. (Q involves all the sites and plaquettes denoted with solid 
circles in this figure. The Ht term contains flavor dependent 
boson hopping between sites and adjacent plaquettes (the Hay 
term in Eq. @), and also between nearest neighbor plaquettes 
(the _ffpp term). (6), the ring exchange term — ti cos(Szz) in 
Eq. pop , which includes eight boson hoppings. 

it separates the Hilbert space into a high-energy subspace 
and a low-energy subspace. The low-energy subspace is 
subject to a local constraint: The sum in the bracket 
in iJ„ vanishes for every Hnk. If we define a symmetric 
tensor £,y as fi^^p = -(-l)'^2(nij^f^ - n), E^.,p^i_^i = 

{—\y{n^^ P^i^l — n) (« 7^ j) this local constraint can be 
compactly written as 

Y.^^S^J=0, (5) 

i 

Here and throughout the paper, Vi denotes the lattice 
derivative: Vi/(r) = f{r+i) — f{r)- Eq. ([5]) simply states 
that £ij is covariantly conserved, or divergence-free. 

The canonical conjugate variables of riij on the lat- 
tice are the phase angles 0ij of the boson creation opera- 
tors: bij ^ e~*^'^ . They satisfy the commutation relation 
[nab, dcd] = iSacSbd, a <b, c< d. Using 6*^ , we introduce 
a symmetric tensor field Aij as An p = —{—lyOa^r and 
Aj.r+i+i = foi^ i j- Under the dis- 

crete symmetries of time reversal T, lattice translations 
T^, along k, and spatial reflection transformations P^, ^, 



Aij — >■ Aij, i ^ k, j ^ k, 
Ai, A,i. (6) 
Notice that we define Pf to be the site-centered reflec- 

k,s 

tion of the lattice, where the origin is located at one of 
the sites of the cubic lattice. The transformation of £ij 
is almost the same as Aij, except that £ij is even under 
time reversal. 

The local constraint on £ij Eq. ^ can be interpreted 
as a Gauss constraint, associated with a partially fixed 
gauge symmetry. It will generate the following gauge 
symmetries on Aij : 

A,j ^ A,j+\7Jj+VjU (7) 

Of course, this gauge symmetry was absent in the micro- 
scopic Hamiltonian Eq. (jj]): It only emerges in the low- 
energy Hilbert subspace with constraints imposed by Hy . 
This mechanism is analogous to the way in which the 
emergent U(l) gauge symmetry appears in the photon 
phase of the 3d quantum dimer models [1, |3| and similar 
spin models 0, Q . If wc interpret Aij as small fluctua- 
tions of the metric tensor gij around the flat background, 

gij ~ % + A^j, (8) 

the emergent gauge symmetry in Eq. ([7]) corresponds to 
the linearized form of the foliation-preserving diffeomor- 
phisms Eq. ([1]) of Lifshitz gravity, with w fi. 

The low-energy dynamics of this system has to be in- 
variant under the gauge transformation Eq. ([7]). The 
lowest order gauge invariant dynamics is generated at the 
eighth order perburbation of t/V, which can be written 
as 

Hofi = ^ — ^ tl COs(i?ijij)p — ^ t2C0s{Rijik)p, 
r i^j i^j^j^k,k^i 

^ijij — ^ —^ijk^kab^kcd^ cAbd, i 7^ j, 
a,b,c,d,k 

Rijik = ^ ejafcefccdVaVcAd, i j,j k,k ^ i. (9) 

With our identification of Aij as the metric fluctua- 
tions around the flat metric, Rijij and Rijik represent 
the six independent components of the linearized curva- 
ture tensor of gij. It is convenient to introduce another 
rank-two symmetric tensor Bij as Ba = ■^eijkRjkjk and 
Bij = —Rikjk- Bij defined this way is also covariantly 
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constant : ViSi. 



0. The transformation of Ba under 



the discrete lattice symmetries of Eq. ^ is the same as 
for Aij. The definition of Bij on the lattice is given in 
the appendix. In geometric terms, Bij are simply the 
components of the linearized Einstein tensor Rij — ^Rgtj 
of the metric 

Now the full low-energy Hamiltonian reads 

r ii i<j 

-^ti cos{Bu)f - ^ t2 cos(Sy )f . (10) 

i i^j 

ti and t2 are both at the order of ^ -P /V'^ . Physically 
ti and t2 terms correspond to high order hoppings of 
bosons. For instance ti term stands for the eighth order 
hopping process depicted in Ref. [13, El and and Fig.©. 
We want to emphasize that the Hamiltonian Eq. (|4]) is 
not fine-tuned, in the sense that one is allowed to turn on 
small local perturbations of any kind that are compatible 
with the global symmetry of Eq. Q , and the low-energy 
Hamiltonian will always take the same form as Eq. (ITUl) . 

The Hamiltonian in Eq. (1101) has a continuum limit 
Gaussian field theory, which characterizes an ABL phase. 
In the field theory we replace the lattice derivative V; 
by the continuum limit derivative di, and expand the 
cosines in iJoS at its minima -~cos{Bij) ^ ^fjl'^ to the 
leading nontrivial order (a spin- wave expansion). After 
we replace Eij by Eij ~ Aij, the resulting Gaussian field 
theory is described by the following Lagrangian, 



2U2 



F F - 
I, while 



The dispersion of the collective modes can be solved 
straightforwardly, and the result is quadratic, implying 
z — 2. Importantly, the linear dispersion is ruled out by 
the gauge symmetry and the lattice symmetries; no other 
terms more relevant than those in Eq. ([TT|) are allowed. 
The terms ^a^jj ^^"^ '^iT^j ^i'-^n principle 

allowed by the symmetry of the system, but it will not 
change the z = 2 dispersion. For instance, ^ 
term can be induced with the V' term in Eq. (| 
X^i^j ^ii^jj *^an be induced with even higher order boson 
ring exchange. 

It is now easy to see that the theory in Eq. ([TT|). to- 
gether with the constraint of Eq. ([5]) , is equivalent to the 
Gaussian limit of the z = 2 Lifshitz quantum gravity 
proposed in [H, Q. Indeed, setting iV^ = as our gauge 
choice in Lifshitz gravity yields Eq. ([S|) as the equation 
of motion, from varying iV^ in the action given in Eq. ^ . 
Furthermore, the most general potential term Sy in the 
z = 2 Lifshitz gravity Lagrangian in 3 + 1 dimensions is 
a sum of two terms, 



i,3 



(12) 



while the kinetic term Sk takes the form in Eq. ([2]). One 
difference between the lattice system and the field theory 
Eq. (fT2l) is that, the lattice system does not have the 0(3) 
spatial rotation symmetry, therefore the low-energy field 
theory deduced from our lattice model does not automat- 
ically acquire this 0(3) symmetry. However, one can tune 
the parameters in Eq. Q to restore the 0(3) symmetry 
in the low-energy limit, and then the field theory will be 
identical to Eq. ([2]) plus Sy- Interestingly, if we tune the 
microscopic parameters to achieve the 0(3) symmetry, 
the ABL phase picks out the special case of the Lifshitz 
gravity models, satisfying the additional property of de- 
tailed balance in the sense of [H, 0]; with a fixed value 
of the coupling A in Eq. ([2]). Extending the relationship 
beyond this simplest case requires that the additional 
terms ^a^jj ^a^jj mentioned above are 

also generated. 

Having established the map bewteen the low-energy 
continuum Hamiltonian of the lattice boson model and 
the Hamiltonian of z = 2 Lifshitz gravity, a few com- 
ments are in order: 

(i) In Lifshitz gravity with the full foliation-preserving 
diffeomorphism invariance of Eq. ([1]), one additional 
global constraint follows from the variation of N{t) in 
the action Eq. ([2|). Since N(t) is only a function of time, 
this yields an integral constraint, equivalent to the van- 
ishing of the total Hamiltonian on physical states. Not 
imposing the 5t — Q[t) invariance as a gauge symmetry 
would eliminate this Hamiltonian constraint. In the lat- 
tice approach, the integral Hamiltonian constraint does 
not seem to be necessary for self-consistency; however, 
we could consider imposing it in addition to the local 
constraint given in Eq. dS]), thus reproducing the full set 
of symmetries of the minimal version of Lifshitz gravity. 

(ii) In [H, 0] , a fully interacting nonlinear version of Lif- 
shitz gravity has been proposed, and it was argued that 
models with z > 1 naturally flow in the IR to z = 1, 
under the influence of relevant terms in 5*^. The most 
relevant terms are R - the Einstein-Hilbert term respon- 
sible for z = 1, and A ~ the cosmological constant term. 
In the self-interacting Lifshitz gravity theory, such terms 
are always expected to be generated under RG, and one 
might wonder why they do not automatically arise in the 
lattice model. The resolution of this puzzle is simple: The 
microscopic structure of our lattice model implies that at 
the z = 2 Gaussian fixed point, the discrete symmetries 
of Eq. ([6]) will hold. It turns out that these symmetries 
not only prevent the relevant terms R and A from being 
generated, they are also incompatible with turning on 
the self-interaction coupling of the full nonlinear Lifshitz 
gravity. 

(iii) One should be more careful with the naive expan- 
sion of the cosine functions appearing in Eq. (fTO|). The 
way we constructed the low-energy Hamiltonian from 
the microscopic degrees of freedom, the Aij variables are 
compact, with radius 27r. This compactification does not 
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have a direct analogy in Lifshitz gravity. Moreover, it 
can lead to topological excitations, which in turn cause 
tunnelling between different minima of the cosine func- 
tions in Eq. Q . Just like in the case of compact QED in 
2 + 1 dimensions, the compactification has the potential 
to destroy the gapless excitations when it is relevant. 

The relevance or irrelevance of the compactification 
and topological excitations can be most naturally studied 
in the dual picture. The low-energy Hamiltonian of this 
ABL phase can be schematically written as H = £^ + B^, 
£ and B are both symmetric, covariantly constant ten- 
sors with six independent components, which suggests a 
self-dual structure exchanging £ and B. The dual field 
hij and the dual momentum tt^ are defined as 

£ij = ^lab^jccN chbd, Bii = '^T^ii, ^i] = T^ij{i j)(13) 

Notice that all the derivatives are lattice derivatives, and 
hij and Tr.y are defined on the sites and plaquettes of 
the cubic lattice (for details, see the appendix). One can 
check the commutator and verify that hij and Tr.y are 
canonical conjugate variables, [7ry,/ifc/] = iSikSji i < j, 
k < I, and we can replace Tr^ by hij. Under this duality 
transformation, £ij and Bij are exchanged, and the dual 
low-energy continuum limit field theory is precisely the 
same as the original model. Therefore, this ABL phase is 
self-dual with the dual gauge symmetry hij hij+difj + 
djfi in the continuum limit. Due to the compactness of 
Aij, £ij and hij both take discrete eigenvalues, therefore 
at the microscopic level the dual Lagrangian allows for 
the periodic potential ^ cos{2TThij) which we refer 
to as the vertex operator. However, this vertex operator 
violates the dual continuum limit gauge symmetry of the 
ABL phase, and hence the correlation function between 
Ov vanishes at distance larger than V/i. Thus, Oy is 
irrelevant at the ABL gaussian fixed point. 

The degrees of freedom in the z = 2 ABL phase indeed 
match the count of degrees of freedom in z = 2 Lifshitz 
gravity [l[ : In addition to the two transverse-traceless po- 
larizations of the graviton, there is an additional "scalar 
graviton" mode. In the lattice model, these three polar- 
izations appear as independent collective modes of the 
Hamiltonian in Eq. ([T0|) . The scalar mode corresponds 
to the scale factor of the spatial metric. It was shown in 
that in classical Lifshitz gravity, the scalar graviton 
can be eliminated by extending the gauge invariance to 
include an anisotropic Weyl symmetry, introduced first 
in [H (and further studied in [l^)- In 3 -f 1 dimensions, 
this requires the dynamical exponent to be z = 3, freezes 
the couphng constant A in the action ^ to he X — 1/3, 
and requires that the spatial part Ly of the action be 
conformally invariant. 

At the microscopic level of our lattice model, the elim- 
ination of the scalar graviton can be arranged simply by 
turning on in Eq. ([5]), When both Hy and Hy/ are 
dominant, the constraint on £ij in the low-energy sub- 



space becomes 

^V,£:,,=0, ^£., = 0. (14) 

i i 

Thus, £ij is now not only symmetric and covariantly con- 
stant, but also traceless. This tracelessness constraint 
can be interpreted as the Gauss constraint associated 
with gauge fixing of another gauge symmetry, which acts 
on Aij via 

Aij -> Aij + S^jif, (15) 

where if is an arbitrary scalar field. Now the trace of Aij 
becomes an unphysical gauge degree of freedom, so there 
are only two gapless modes at low energy. 

The new gauge symmetry in Eq. (jisp will turn out 
to be the anisotropic Weyl invariance of Lifshitz gravity 
mentioned above, here in the linearized approximation 
around the Gaussian z = 3 fixed point. In order to see 
that the new constraint is forcing the model to z = 3 at 
low energies, note that the low energy Hamiltonian is also 
modified by turning on V' . Let us assume the dynamical 
term can be written as 

H^s' = X! X! cos{aCii)f - ^ t4 cos{bCij)f, (16) 

r i i^j 

where Cij{Aki) is a linear functional of A^i. This new 
tensor Cij must be invariant under the gauge transfor- 
mations of Eq. ([7]) as well as Eq. p^ . The tensor of the 
lowest order in derivatives satisfying this requirement is 
the (linearized form of the) Cotton tensor [1113]: 

C^J=e^klyk{RJl-\R5■Jl). (17) 

Here Rji = J^k^jkik, R = J2j Rjj the linearized 
Ricci curvature and scalar curvature respectively. Un- 
der the lattice symmetries of Eq. ([6]) , the Cotton tensor 
transforms as 

T . Cij y ^ijj 

Cij — >■ —Cij, i ^ k,j ^ k, 

Cii —Cii. (18) 

One can straightforwardly verify that Cij is gauge invari- 
ant, symmetric, traceless and covariantly constant: 

C^j = J2 C« = 0' II ^^^-J = 0- (19) 

i i 

The definition of Cy on the lattice is given in the ap- 
pendix. From our microscopic Hamiltonian, the effective 
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low-energy term -ffoff' can only emerge at fairly high or- 
der of perturbation. For instance, the ^3 term can be 
generated at order 32 in i/V . Now the full low-energy 
Hamiltonian reads 

r ii i<j 

- ^ ^3 cos{Cii)p - ^ ^4 cos(2Cij)f . (20) 

i i^j 

Since Cij involves three spatial derivatives, after the spin- 
wave expansion — cos{bCij) ^ the ABL phase 
which is described by the continuum limit Gaussian field 
theory of Hamiltonian in Eq. (|20p has collective excita- 
tions with the z = 3 dispersion. Notice that the lattice 
symmetry does not require ^3 = t^, but the ratio ts/t^ 
can be tuned by the ratio ^1/^2 on the lattice. 

In 0] , a nonlinear self-interacing z = 3 Lifshitz gravity 
has been constructed. In this construction, the nonlinear 
Cotton tensor plays a central role, with Sv oc CijC^^ . 
When ^3 — t4, the low-energy Hamiltonian of the ABL 
phase is identical to the Gaussian point of the z — 3 
Lifshitz gravity of This emergent theory is gauge 
invariant not only under linearized foliation-preserving 
diffeomorphisms, but also under the z = 3 version of the 
anisotropic Weyl transformation [l|, 0, [3| • 

Are the gapless z = 3 modes of this model stable 
when Aij are treated as compact variables? To address 
this question, it is again convenient to go to the dual 
formalism. The continuum limit Hamiltonian takes the 
schematic form H = £^ + C^, and the fact that £ij and 
Cij are both symmetric, traceless and covariant tensors 
strongly suggests this theory also has a self-dual struc- 
ture. To prove the self-duality, we define dual variables 
hij and as 

£ij =Cij(hki), TTu = -Cu{Aki), TTij ^C.ij{Aki), i ^ .(21) 

Here we have treated Cij as a linear functional of the 
tensor field Aij or hij. Unlike hij defined previously, 
ha are defined on the dual lattice sites, which are the 
centers of the unit cubes, and hij are defined on the dual 
plaquettes, which are the links of the original lattice (for 
details, see the appendix). For instance, a link along the 
z direction is a dual plaquette in the xy plane. Again, 
one can straightforwardly verify that hij and liij are a 
pair of conjugate variables, and £ij and Cij are exchanged 
under duality. The dual graviton hij of this ABL phase 
enjoys the same gauge symmetry as Aij, hence the vertex 
operator ~ cos{2'nhij) is irrelevant. Therefore the 
2 = 3 ABL phase is also a stable gapless phase. 

This self-duality structure completes the argument of 
the stability of both the z = 2 and z — 3 ABL phases 
in this model. The self-duality can also be proved in the 
Euclidean space-time, in the same way as the duality of 
ordinary classical statistical mechanics models . If 



the dual theory of a lattice gauge model did not have a 
large enough gauge symmetry, one would have to fine- 
tune the system to get a gapless ABL state. The most 
well-understood example is the compact QED in 2 -|- 1 
dimensions, where the dual theory is a U(l) rotor model 
without gauge symmetry. In that case, the vertex op- 
erator which corresponds to the monopoles in spacetime 
will destroy the gapless photon phase. Another example 
studied recently is the quantum plaquette model with 
gauge symmetry Aij — >■ ViVj(p. The dual theory of the 
plaquette model does have a gauge symmetry, but this 
symmetry is not strong enough to protect the gapless 
ABL phase |:20:]. 

In Ref. |21j], graviton- like collective modes were also 
obtained through a quantum boson model on the lattice, 
but the relation of such lattice models to Lifshitz gravity 
was not noticed. It is worth emphasizing that in our con- 
text, the Lifshitz-type graviton Hamiltonian in Eq. ((20)) 
can be derived without any fine tuning from a simple bo- 
son model of Eq. Q through perturbation theory, and 
both the z — 2 and z = 3 ABL phase are stable due to 
their self-dual nature. 

Besides the difference in the dispersion, the z = 2 and 
z = 3 phases also have different algebraic correlations. 
For instance, in the z — 2 phase the equal-time density 
fluctuation correlation falls off as 

{Sn{0)dn{f)) ^ {B{h)aB{h)p) ^ 1, (22) 

while for the z — 3 phase this correlation falls off as 

{6n{Q)6n{r)) ^ {C{h)oC{h)p) ^ 1. (23) 

(In these two equations, the flavor dependence of the cor- 
relation has been ignored for simplicity.) 

In the original Hamiltonian, the z = 2 and z — 3 
phases were obtained by dialing a small versus large value 
of the coupling V' in Hyi . As a result, by increasing 
V' one can drive a phase transition between these two 
ABL phases with different dynamical scalings. Com- 
pared to the z = 2 phase, the z = 3 phase has one ex- 
tra gauge symmetry, of anisotropic Weyl transformations 
P^ . Therefore, if we start with the z = 3 phase, this 
phase transition will be reminiscent of the Higgs transi- 
tion. In the quantum dimer model, the U(l) gauge sym- 
metry is Higgsed by condensing the dimer vacancies that 
carry U(l) gauge charge. The phase transition can be de- 
scribed by the rotor Hamiltonian H,. = —t cos( V6' — A) 
[2^ , where 9 is the phase angle of the dimer vacancy cre- 
ation operator tp ^ e'^^ . After the condensation of tp, the 
gauge field A^ acquires a longitudinal mode by absorbing 
the Goldstone mode of '0, which makes the gauge field 
an ordinary gapped vector field. 

In our case, the transition between the two ABL phases 
can also be intuitively described as closing the gap of the 
trace mode of £ij , which can be described by condensing 
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(f in Eq. ([15]). Since all the gauge symmetries in Eq. ([7]) 
are still preserved after the transition, the condensate 
of ip should not violate Eq. ([7]). With these observa- 
tions, this transition can be described by the following 
Lagrangian, 



L = -ip^-^ti cos(V^^ + V^^ - R^j^j) 



t2 cos(V-,- Vfcvs - Rijik) 



(24) 



The condensation of 95 changes the spectrum from two 
z = 3 modes of Aij and one gapped mode of Lp to three 
z = 2 modes of Aij , in a generalization of the ordinary 
Higgs transition to the scalar mode of Lifshitz gravity. 
When if is ordered, this Lagrangian restores the ring ex- 
change terms in Eq. ([TU]) of the z = 2 phase; when Lp 
is disordered, after integrating out one should recover 
the phase which is invariant under the transformations in 
Eq. (O and Eq. The z = 3 ABL phase is then the 

only candidate. From the condensed matter perspective, 
the phase transition in Eq. (j24p is beyond the ordinary 
Ginzburg-Landau paradigm, because neither one of the 
two phases can be characterized by a local order parame- 
ter. More thorough RG studies for Eq. (|24|) are required 
to determine the nature of this transition. 

One interesting question to ask is whether we can ob- 
tain relativistic gravitons, with a linear dispersion, from 
the lattice. In Ref. [llj , it was proposed that a long-range 
interaction can change the dispersion to z = 1, but a lo- 
cal theory leading to z = 1 gravitons at long distances 
is still unavailable. As was discussed in Ref. ll| and 
Ref. [2^, a Chern-Simons like term AijBij can lead to a 
linear dispersion, but this term is only gauge invariant up 
to a boundary term; therefore, it cannot be generated in 
the same way as the ti and ^2 terms in Eq. (fTO| through 
perturbation theory. One possibility is to generate this 
term by coupling the graviton field Aij to a matter field 
with a gapless boundary state, just like the CS term for a 
U(l) gauge field can be generated by couphng the gauge 
field to a massive Dirac fermion with edge states. 

In addition to the linear dispersion, another meaning- 
ful goal is to obtain the full nonlinear Lifshitz gravity of 
self-interacting gravitons, instead of the linearized the- 
ory at the Gaussian fixed point. As we pointed out in 
our comment (ii) below Eq. (|12p . this difficulty is inti- 
mately related to the existence of the discrete symmetries 
in Eq. ([6]), implied by the microscopic dynamics of the 
lattice model. These symmetries do not allow the nat- 
ural self-interaction coupling of Lifshitz gravity P, 0] to 
be turned on. It is an interesting challenge to see if our 
framework can be extended so that its discrete symme- 
tries no longer prevent the self-interaction of gravitons. 
Note that the discrete symmetries of Eq. ^ act natu- 
rally on the Aij variables representing the fluctuations 
around the fixed flat background, but they do not ap- 



pear to have a geometrically natural extension to the 
full metric gij. Thus, they are associated with the fixed 
flat spatial geometry, here represented by the fixed flat 
fee lattice. These background-dependent discrete sym- 
metries indeed played an important role in our construc- 
tion of Lifshitz gravity from a lattice system: They pre- 
vented the Einstein-Hilbert term and the cosmological 
constant term from being generated, allowing the z = 2 
and z = 3 ABL phases to be stable at low energies. Since 
the full nonlinear Lifshitz gravity of [ll[3| does not require 
a choice of a preferred flat background, it is natural to 
speculate that attempts to turn on the self-interaction of 
gravitons in our lattice framework may ask for the un- 
derlying lattice itself to become dynamical. We will leave 
these topics to future studies. 
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APPENDIX: TENSOR FIELDS ON THE LATTICE 

On the lattice, tTzz on site f is defined as 

'^'^zz^f ( 1) {^yy ,r-\-x ^" '^^yy.r^ ^yy.r—x 



^xx^r—y 



2^ — , X I il 20 — X , y 

xy,r+^ + ^ xy,r-^ + ^ 



-2^xy,r+f-i '^Qxy,r^i-^. 



Bzz- (25) 



The ring exchange coa{Bzz) corresponds to the high order 
of boson hopping depicted in Fig. [2). iTyz on plaquette 
r+ ^ + ^ reads 



yz,r+^ + ^ 







( ^yz,r+x+^ + ^ 




~t~ ^xx^r ~^ ^xx,f- 




^xy,f+§ + ^ 




^ xy,r+z+^ + ^ 


xy,r+z—^ + 






^xz,r+y+^ + ^ 


~ ^xz,f+y-^ + 



+ 1 



-Byz. (26) 

The tensor field defined this way satisfies J2i '^iBij = 0. 
The relation between £ij and hij is the same as that 
between Bij and Aij . 
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The dual variable tth is defined on the dual sites, which 
are the centers of the cubes, and -iry with i ^ j is located 
on the links r + |, with k ^ i, k ^ j. Wxx on the dual 
site r+l + l + l is defined as 

+^a;y,f+z+f + | ~^xj/,f+f + | =^^^- (2'') 



■Kxy on the link r + | is defined as 



"I" ^zz,r ^xx,f+z "I" ^xx,f — 2Ca;y. (28) 

Cij defined this way is symmetric, traceless, and covari- 
antly constant. The relation between £ij and the dual 
variable hij is identical to that between TTj^ and Aij , after 
exchanging sites with cubes, and plaquettes with links. 
On the lattice, the divergence of Cij reads 



C ' ■ — C ■ - ■ 

a;x,f+^ + -| + | xx,r—^ + ^ + -^ 



A- C - — C 

~ ^xy,f+y+^ *^a;j/,r+f 



+ e ^^.^y -c 



(29) 
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